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SUMMARY 


Fatigue crack growth associated with cyclic (secondary) plastic flow near a crack front is 
modelled using an incremental formulation. A new description of threshold behaviour under small 
load cycles is included. Quasi-static crack extension under high load excursions is described using 
an incremental formulation of the R-(crack growth resistance)- curve concept. 

The integration of the equations is discussed. For constant amplitude load cycles the results will 
be compared with existing crack growth laws. It will be shown that the model also properly 
describes interaction effects of fatigue crack growth and quasi-static crack extension. 

To evaluate the more general applicability the model is included in the NASGRO computer code 
for damage tolerance analysis. For this purpose the NASGRO programme was provided with the 
CORPUS and the STRIP- YIELD models for computation of the crack opening load levels. 

The implementation is discussed and recent results of the verification are presented. 


INTRODUCTION 


For over two decades models of fatigue crack growth have been based on empirical laws that 
relate the amount of crack growth in a load cycle to the stress intensity factor range AK = K max - 
K^ m i n or the effective [1] range AK eff = - K . Correction factors were included for near 

threshold behaviour and accelerated growth in the nigh K regime. 

From a physical point of view such crack growth laws are speculative because crack growth and 
plastic deformation are irreversible processes that depend on the loading history. By nature, such 
processes must be described in an incremental way and properly integrated to obtain the amount of 
crack growth for a load cycle or the part of a load cycle for which the incremental description is 
valid [2, 3]. Clearly, such a new description allows that a distinction is made between the part of a 
load range where secondary (cyclic) plastic flow is observed and the part where primary plastic 
flow developes under monotonic increasing loads. For each of these domains an incremental crack 
growth law can be formulated. Then after integration over the appropriate load ranges the 
contributions to the crack growth rate for the load cycle under consideration are obtained. In a 
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similar way "range pair" (or "rain flow") principles may be used to select the appropriate crack 
growth equations and associated ranges of applicability. In addition, the incremental formulation 
allows the introduction of other terms representing time and/or environment dependent crack 
growth. In these applications the integration of the incremental equations becomes more 
complicated, however, with the increased capabilities of both the numerical techniques and the new 
generation of computer systems numerical integration appears to be feasible. 

In this paper the formulation of two incremental crack growth laws, one describing fatigue crack 
growth and a second one for static crack extension, is discussed. For constant amplitude loading the 
equations are integrated analytically to obtain the crack growth rate per load cycle. The result can 
be compared with the common crack growth laws and allows the material parameters to be 
determined from the results of simple fatigue crack growth tests. For more complicated load 
sequences like block programme loading and arbitrary, cycle by cycle defined sequences, a method 
for integration of the crack growth rate is discussed. The models and methods were implemented 
into the NASGRO (ESACRACK) software [4] and thoroughly tested. In view of uncertainties about 
the values of constraint factors a preliminary verification was executed using some new test results 
obtained for three materials: a titanium alloy Ti-6A1-4V, the aluminium alloy 7075--T73 and the 
COLUMBUS" skin material 2219-T851 for the European space station manned module. Two 
material thicknesses were involved: 2 or 4 mm sheet and, 10 mm, plate material. The initial 
precracks were edge cracks, surface cracks and through the thickness cracks. After preliminary 
verification the models were readjusted and the constraint factors specified in more detail. Some of 
the proposed modifications are discussed. The final verification also includes a large number of test 
cases collected from the open literature in addition to in-house test results [5]. 

This paper describes the status of the study after completion of the preliminary verification and the 
model evaluation. 


THE LOADING AND CRACK GROWTH REGIMES DISTINGUISHED IN THE UPWARD 

PART OF A LOAD CYCLE 


In the upward part of a load cycle different regimes can be distinguished, depending on the 
characteristics of the plastic deformation behaviour and the state of opening of the crack. To 
illustrate these domains in figure 1 the loading path is shown in a K versus c plot, where c is the 
crack length. The different loading regimes are indicated and discussed one after another. 

Closed crack regime 1, K min < K < K 

Starting at the minimum stress intensity factor the load is increased until the crack opening 
level K op is obtained. In this first regime 1, characterized by <_ K < K op , the crack is at least 
partly closed and the contact areas on the crack surfaces decrease when the applied load is 
increased. Although the stress intensity factors in this regime are calculated assuming the presence 
of the crack, it is clear that the effective loading of the crack tip region is very small and no crack 
growth is assumed in this regime. 

Opened crack but no growth regime 2, K op <_ K < K op + 5 K th 

At level K op the crack is fully opened, but, it takes another increase by to initiate crack 
growth. Obviously, some crack tip blunting occurs in this regime 2. Models and empirical equations 
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for computation of values for K op and are discussed later on. 


Fatigue crack growth in regime 3, K + bK th < K < K+ 

Upon a further increase of the applied load crack growth is initiated when the stress intensity 
factor K exceeds the level K op + SK^. In this regime 3 alternating secondary plastic flow is 
observed in a relatively small plastic zone. At level K = K*, however, the cyclic plastic flow is 
assumed to change to primary plastic flow. This transition is characterized by a discrete jump in 
plastic zone size and a loss of load history effects on the state of deformation. To describe the crack 
growth behaviour in regime 3, corresponding to K op + < K _<. K*, the following incremental 

crack growth law is adopted 


dc 


C[(K-K op ) n + C 2 5K t £(K-K op ) n -P ]dK 


( 1 ) 


In this expression the first term on the right hand side is an incremental form of Elber’s law. The 
second one is added to describe threshold effects, if present. The power n-p follows from the 
requirement that the units of both terms must be the same. 

At initiation of crack growth, when, K — K op + SK^, it follows that 

[— | = 5KJJ(C 1+ C 2 ) (2) 

i. dK J th 



Fig. 1 Different loading and crack growth regimes in one (half) load cycle 
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Obviously, there exists a relation between the material parameters and C 2 , the threshold level 
and the slope of the crack growth curve (see Fig. 1). At the present time knowledge about this 
slope is lacking. Therefore, a convenient value for the slope is selected to simplify the equations. 
This value is 


dc j 
dK j th 


pC,5K^/(n + l) 


(3) 


To obtain the amount of crack growth Ac s created under regime 3 the crack growth law (1) must be 
integrated over the range K op T < K <_ K*. There results 


Ac 


s 



n + 1 


K op ) 


n + 1 


1 - 


-o 

— 1 


K - K op 



(4) 


In the next section this equation is used to derive an equation for computation of the transition level 
K*. If the transition level is above K max then the upper bound in the integration of eq. (1) is K max 
and in eq. (4) K* is substituted by K max . 

Quasi-static crack extension regime 4, K* < K < K niax 

Loading above the transition level K* is assumed to induce quasi-static crack extension. In this 
regime the plastic deformation behaviour takes place under monotonic increasing loads. This implies 
that the effects of secondary cyclic loading on the actual material behaviour are lost. Thus, the 
crack opening load and threshold behaviour become insignificant [2, 3]. Moreover, the plastic zone 
sizes are much greater. To describe crack growth in this domain we will adopt the incremental 
formulation of the R (or J) curve approach. Assuming small scale plastic behaviour and small 
amounts of static crack extension the crack growth law adopted is written as 

dc = C p K m dK (5) 

In addition, for cases where wide scale plastic deformation occurs or the amount of static crack 
extension becomes large we may choose to introduce new -or sub- regimes and formulate the 
applicable crack growth law in such a way that it describes these processes properly. 

The incremental crack growth law must be integrated over the applicable range to obtain the 
contribution Ac p to the crack growth increment for a load cycle. There results 

Ac„ = -El. [k- 1 - kT 1 ] (6) 

m + 1 L J 


DETERMINATION OF THE TRANSITION LEVELS K^, K op + SK^, AND, K* 


The crack opening level , K op 

In this study the following crack opening models are used: The mechanical STRIP- YIELD (= 
modified discretized Barenblatt/Dugdale) model [6, 7, 8, 9] and the empirical CORPUS 
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(Computation Qf Retarded crack propagation Linder Spectrum loading) model [10, 1 1, 12]. Both 
models were included in the NASGRO software [4, 5]. In both crack opening models constraint 
factors are introduced to account for 3D effects on the yield stresses of material in front of the 
crack tip (in tension and in compression) and the yield stress for loading of material in compression 
in the wake of the crack. These effects are thought to be determined by the state of stress, the 
uniaxial yield limit a and the sheet or plate thickness T. Only two parameters are considered to be 
material parameters, namely: the ratio a NEW of the cyclic yield limit in tension over its value in 
compression, and, the uniaxial yield limit a. The modelling of the constraint effects will be 
discussed later on. The STRIP- YIELD model included in NASGRO was derived from the NLR 
model [7]. Newman’s method [13] for computation of the crack opening load from the contact 
stress solution at the minimum load in the cycle was adopted. 

Since in some of the applications the static crack extension plays an important role its effect on the 
crack opening level K must be accounted for. In the STRIP-YIELD model this is simulated by 
unzipping elements in front of the crack tip one after another. In this study unzipping takes place at 
the minimum load in the cycle (no additional plastic deformation is added to the wake). In the next 
version of the model, currently under development, the elements are unzipped at the proper stress 
intensity level. Then, the effect of the additional stretch resulting from static growth of the crack is 
automatically included. 

In an application of the CORPUS model [10] the crack opening level is calculated using the 
crack opening function f(R, c* NEW , S inax /o). This function is chosen in such a way that the 
STRIP- YIELD analyses results of constant amplitude data are described accurately. The CORPUS 
model defines a set of rules for selection of the representative load ratio R and the stress level S max 
from the loading history. In the next version of CORPUS implemented in the NASGRO software the 
crack opening function will be corrected for the effects of static growth and other improvements of 
the STRIP- YIELD model. Then, the model can be used to determine the material parameters also 
for the STRIP- YIELD model in an efficient way. 

To illustrate such corrections on the crack opening function the correction accounting for the 
effect of static crack extension is discussed in a similar way as in reference [8] where a first order 
approximation was used. Thus, the crack opening level is written as 

K op / K ra ax “ f(R> a NEW>S max / °) ^ 


From the result of eq. (7) we can calculate the crack opening load level S op and estimate the 
fictitious, STRIP- YIELD based, contact area Ac w near the crack tip at the minimum load level 


Ac w 


c 


f 

1 - sin 

! _ 



f s ~S • 1 

- 

- 

T 

1 „ °p min 



2 

c 

cr„ 




L y J 


j 


( 8 ) 


where Is the yield stress of the material in the wake of the crack. Assuming that the contact area 
is extended by the amount of growth Ac p it follows for the correction on the crack opening stress 
that 
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AS 


2 <r 


c-Ac. 


op 


arcsin 


in eq. (9) c denotes the crack size. 


- arcsin 


. ] 

- 

c-Ac w 


c-Ac 


P J 

— 


(9) 


The effective threshold level, K + 8K (h 

The effective threshold range is derived from the empirical relation used in the NASGRO 
software. There results 


«** = 


1 NEW> S max/ a ) 


NEW > J max‘ 

1 - R 




( 10 ) 


where - K min represents the NASGRO threshold stress intensity range and f(R,a NEW ,S max /r) is 
the opening function. The threshold value applies to constant amplitude loading of through the 
thickness cracks. In the applications to part through cracks and to variable amplitude loading the 
value of SKfa is reduced. 


The stress intensity factor at transition from cyclic (secondary) plastic flow to primary flow, K m 
To derive an expression for computation of the level K* at which the transition from secondary 
plastic flow to primary plastic flow occurs a loading sequence as given in figure 2 is used. The first 
spike load 1 has created an overload plastic zone and at the time the second spike is applied the 
crack tip is assumed to be still situated inside this overload plastic zone. In terms of K versus crack 
length the situation is sketched in figure 3. Then, from a geometrical consideration and the 
assumption that the transition occurs when the actual primary plastic zone D* touches the end of the 
overload plastic zone D SP the following expression can be written 

c “ Cg P + Ac s + D* - D sp = 0 (II) 


where c is the crack length at initiation of crack growth in the actual overload cycle 2. The 
secondary crack growth increment Ac s , given by eq. (4), is a function of K*. 
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Fig. 2 The loading schedule used to derive an equation for computation 
of the transition level from secondary to primary plastic flow 






stress intensity 
factor K 


C “ C sp + ^ C 8 + “ ^sp 

where: 


* 8a 2 aJ 
, *<K 8p ) 2 

; sp = _ _o o~ 

8 a 2 a|p 


amount of crack 
growth since 
application of 
the first spike 


spike plastic zone 
at K = K cn 


crack size c 


actual plastic 
zone at K = K* 


spike plastic zone size D S p 


Fig. 3 Geometrical consideration used for determination of the transition level K * 
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Further, the plastic zone sizes will be approximated using a First order estimate, that is 


and 


D* 


x K » 2 
8 (a*<0 2 


D 


SP 


x ( k sp ) 2 
8 Kp ") 2 


(12) 


(13) 


where a* and or sp are constraint factors introduced to account for 3D effects on the yield stress for 
primary plastic flow. In the preliminary verification the values a* and a sp are set equal to the 
material parameter of NEW . 1° next version of the model a more correct description of the 
constraint behaviour will be introduced. From eqs. (11) to (13), and eq. (4), the value of K* can be 
solved in an iterative way. Finally, the value calculated for K* must be compared with K max in 
order to establish the presence of primary plastic flow and associated quasi-static crack extension. 
Primary flow is absent, if K* > K max . 


APPLICATION OF THE INCREMENTAL CRACK GROWTH LAWS 


Constant amplitude loading 

From the foregoing considerations it is concluded that the crack growth rate per load cycle 
Ac/AN can be written as 


Ac/AN = Ac. + Ac n (14) 

where Ac s and Ac p denote respectively the amounts of crack growth associated with secondary 
plastic flow and with primary plastic flow. After substitution of eqs. (4) and (6) and using the 
knowledge that in constant amplitude loading K min <_ K + 8K ^ < K* _< K we can calculate 
Ac/AN. 


In constant amplitude loading K* is slightly lower than K max in the major part of the crack 
growth curve. Using this result and recognizing that, in figure 3, c sp can be identified with c eq. 
(14) can be approximated by 


Ac/AN 


CflWK,,,,) 1 ’* 1 


■ 


' 





1 - 
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P j 


1 + 

K max ' 

m-2 

^max~^op 

J 


* 

- 


(15) 


where C — Cj/(n+ 1) and the parameter K ref is a combination of the constraint factor a* and a 
number of material parameters according to 


(1/K re f) m_2 = C p mi(a t o) 2 /2. 


(16) 
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In eq. (15) the first part on the right hand side is equal to Elber’s crack growth law. The second 
factor accounts for threshold behaviour in a common way and the third, non-singular, part accounts 
for static crack extension at higher crack growth rates. Eq. (15) is very similar to common crack 
growth laws like the NASGRO modified Forman law. This allows the conversion of parameter 
values by requiring that the different multipliers in eq. (15) have approximately the same effect on 
Ac/AN as the corresponding factors in the crack growth law under consideration. It is concluded 
here that the incremental description of crack growth has given a physical interpretation of the 
increased crack growth rate at high K max levels. Further, it is noted that, for any other initial slope 
of the crack growth curve an expression similar to eq. (15) can be derived. 

Constant amplitude loading interrupted at regular intervals for application of a spike load 
excursion 

Relatively large amounts of static crack extension are observed during spike loading excursions 
applied at regular intervals in an otherwise constant amplitude load sequence, as indicated in figure 
4a. Such sequences can be used to determine the material parameter values involved in static crack 
extension by fitting the integrated crack growth equation to measured data points. In general, an 
overload ratio OL is chosen of the order 1.5 to 1.7. The number of constant amplitude cycles in 
one block is in the order of 500 to 2000. 




Fig. 4b The spike load sequence for the case of non-interacting spike loads. (Large N CA ) 






In such sequences static growth is absent in the constant amplitude cycles provided that 
K max < K* _< OL * K max . Thus, the crack growth equation (4) for secondary (cyclic) plastic flow 
applies to the constant amplitude cycles. For the spike load exclusion K* must be solved iteratively 
from eq. (1 1) and it must be verified that K* _> K max . Then, it can be derived that the average 
crack growth rate resulting from application of one block of N CA constant amplitude cycles plus one 
spike load excursion can be written as 


Ac/AN = 


c [(Kmax-Kop)"* 1 * N CA ♦ (K,-K op ) n *>] ♦ [ K ™ + 1 - K^'] 

(N ca + 1) 


(17) 


where it is assumed that is independent on the crack size and can be calculated using the simple 
crack opening function (7) when S max , a NEW , and R are associated with the spike load cycle. 

In the case K* <. K max the contribution of static growth of the crack must be calculated. This 
situation is shown in figure 4b. To calculate the crack growth rate for this part of the constant 
amplitude cycles eq. (15) can be used. The result can be added to eq. (17) as a weighted average to 
obtain crack growth rate for one block. * c 

In the applications discussed in this paper the Overload Level OL is chosen to be so high that 
threshold behaviour is absent, and, K* > K inax . 


Block programme loading and randomized or cycle by cycle defined load sequences 
In the aerospace industry load spectra are often defined as blocks consisting of a number of 
different load steps each of which containing a distinct number of constant amplitude cycles of a 
given amplitude and mean load level. In the analysis the load steps can be applied one after another 
and, using some additional assumptions about the transition from one load step to the next one, we 
can analyze this specific load sequence derived from the load spectrum (see Fig. 5). Alternatively, 



time ► 

Fig. 5 Three load steps from a block programme load sequence . Parts where crack growth during 

primary plastic flow (PPF) occurs are indicated. Note the different behaviour in the first cycle of 
a load step compared to the remaining part of the constant amplitude cycles. 
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we can randomize the cycles in such a way that, after verification by application of counting 
methods, the sequence properly represents the load spectrum. In this way different sequences can be 
constructed all representing the same load spectrum. It is noted, however, that the load sequences do 
not inevitably represent the same fatigue loading experience. 

Characteristic for the aircraft wing and tail load sequences is the presence of air-ground-air load 
cycles that are included in the load sequence at the proper intervals. In general, these sequences 
contain many (short) periods of constant amplitude loading that can be described as load steps. 

Other cycles can also be described as load steps by defining them as single (half) cycles. 

In the NASGRO software such a definition of load steps was applied. Each first cycle in a load step 
is analysed using a separate crack growth prediction module that calculates the crack growth 
increments Ac p and Ac s for that particular cycle. 

The remaining constant amplitude cycles are analysed by taking discrete steps Ac and (using mid 
point integration at c + l AAc) the number of cycles corresponding to the stepsize Ac is calculated 
and subtracted from the number of cycles left in the load-step. In this way each of the load steps is 
analysed one after another. 


PRELIMINARY VERIFICATION AND MODEL EVALUATION 


Model definition 

In an application of the crack growth and crack opening models additional assumptions must be 
made to quantify constraint effects on the yielding and crack opening behaviour. 

In an early stage of this study it was decided to execute the preliminary predictions using a highly 
simplified system of constraint factors. The yield stresses adopted are illustrated in figure 6. The 
material yield parameter cx NEW is accounting for different yielding behaviour of material loaded in 
tension compared to material yielding in compression. In compression the yield stress for the wake 
of the crack is assumed to be the same as the yield stress of material in front of the crack tip. All 



Fig. 6 Yield stresses used for the preliminary predictions 
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other constraint or Bauschinger effects are assumed to be absent. Later on, after completion and 
evaluation of the preliminary predictions a more correct definition of constraint factors will be 
introduced. In the conclusions and discussion some directives are given for such a system. In the 
application of the CORPUS model the NLR crack opening function is used. The function is given in 
appendix A. 

Test programme 

An experimental programme was defined and executed to collect crack growth data for 
determination of the material parameters in the crack growth law and to verify the models as 
implemented in the NASGRO software. 

The materials involved are Ti-6A1-4V, Al-7075-T73(5 1) and A1-2219-T851. Both sheet and plate 
materials are used. The sheet specimens were centrally precracked (through the thickness crack 
2c 0 = 6 mm). The specimen width is 120 mm in all cases. The 10 mm thick plate specimens 
contained a comer crack or a surface crack. The crack size was 1 mm (a^Co = 1). The load 
sequences applied in the test are constant amplitude loading at R = 0.05 and R = 0.70, repeated 
and single spike load sequences and a Space Transportation System STS (pay load) spectrum. The 
STS spectrum is applied in blocks of constant amplitude loads (full cycles). For the Al-7075- 
T73(51) material the test programme is given in tables la and lb. For the other two materials a 
similar test program was executed. As the conclusions drawn from the results of the preliminary 
verification are the same for all three materials, in this papier the discussions are concentrated on the 
7075-T73(51) material. In reference [5] the results obtained for the other materials are discussed in 
detail. 

Table la Test programme for 2 mm thick 7075-T73 sheet material. Specimen width is 120 mm 


Group 

Spec. 

Fatigue loading programme 

Basic loading 

Spikes 

I 

SO 17075 
S027075 

C.A., R = 0.05, S,_ v = 89.2 MPa 
C.A., R = 0.05, S milv = 89.2 MPa 

: 

n 

S037075 

S047075 

C.A., R = 0.7, S max = 282.5 MPa 
C.A., R = 0.7, S fflax = 282.5 MPa 

• \ 

in 

S057075 

S067075 

S077075 

S087075 

C.A., R = 0.05, S max = 89.2 MPa 
C.A., R = 0.05, S max = 89.2 MPa 
C.A., R = 0.05, S max = 89.2 MPa 
C.A., R = 0.05, S max = 89.2 MPa 

S sp — 151,6 MPa, AN = 500 cycles 

S = 151.6 MPa, AN = 1000 cycles 

S = 151.6 MPa, AN = 2000 cycles 

S — 151.6 MPa at c = 5, 7, 10 mm j 

“P 

IV - 

SI 17075 
S 147075 

C.A., R = 0.05, S mQV = 101.1 MPa 
C.A., R = 0.05, S mov = 101.1 MPa 

111 HA 

S = 151.6 MPa at c =5, 7, 10 mm j 

*P 1 

S g p = 151.6 MPa, AN = 2000 cycles 

v 

SI 27075 
SI 37075 

STS spectrum, S max = 151.6 MPa 
STS spectrum, S max = 282.1 MPa 

: 
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Table lb Test programme for 10 mm thick 7075-T735 1 sheet material. Specimen width is 50 mm 


Group 

Spec. 

Type of 
notch 

Fatigue loading programme j 

Basic loading 

Marker loading or spikes | 


PO 17075 

■ 

— 

■ 

C.A., R = 0.05, = 135 MPa 

Marker loading: C.A., R = 0.6 

S max " 66 8 M P a at c = 5, 7.5, 10, 13 mm 


nu 

SC 

SC 

CC 

C.A., R = 0.05, S max = 135 MPa 
C.A., R = 0.05, S max = 135 MPa 
C.A., R = 0.05, S max = 135 MPa 

S = 229.5 MPa, AN = 2000 cycles 
S = 229.5 MPa at c — 5. 7, 10 mm 
S S p = 229.5 MPa, AN = 2000 cycles 

1 

H 

■ 

■ 

C.A., R = 0.05, S max = 153 MPa 
C.A., R = 0.05, S ?viav = 153 MPa 

indA 

S sp = 229.5 MPa at c — 5, 7, 10 mm j 

S sp — 229.5 MPa, AN = 2000 cycles 1 

n 

P067075 J 

sc 

STS spectrum, S — 340 MPa 

IU<IA 

1 


c = half crack length along specimen surface 
SC — Surface Crack 
CC = Comer Crack 


Determination of the material parameters 

The material parameters appearing in the crack growth laws (4) and (15) can be determined by 
minimizing the distances of the measured data points to the plotted graph of the crack growth 
equation in the log AK = K max - K min versus log (da/dN) domain. From the test matrix defined for 
centrally cracked sheet specimens the following tests were chosen for determination of the material 
parameters: 

SO 17075 Constant amplitude loading at R = 0.05. 

S037075 Constant amplitude loading at R = 0.70. 

S067075 Constant amplitude loading, interrupted each 1000 cycles for application of a spike load. 
The spike load ratio S sp /S max = 1.7. The load ratio R = 0.03 (for the spike). 

The test results obtained for these 3 tests are collected in one series of dc/dN(i) - AK(i) data 
sets. Then, using a standard routine the parameter values in the crack growth law are determined in 
an iterative way such that 

sum = Y, [ ln(dc/dN(i)) - In F(AK(i), K op (i), S max> R) ] 2 (18) 


is a minimum. Thus, the least squares fit is applied on a log scale. In eq. (18) function F(AK(i), 
K 0 „(i), S max , R) represents the right hand side of the crack growth law eq. (15) for S017075 and 
S 037075. In the evaluation of test S067075, eq. (17) is used to calculated the average crack growth 
rate. 
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The crack opening loads K op (i) and, therefore, the values of the material parameters in F depend 
on the crack opening model that is chosen. In the determination of the material parameters the 
CORPUS model is applied in combination with the NLR opening function. A complicated iterative 
scheme based on a STRIP- YIELD analysis of K op (i) was thought to be impractical. The "fit" results 
are indicated in figure 7b. For comparison the fit results obtained for Ti-6A1-4V and AI-2219-T851 
are given in figures 7a and 7c. From figure 7 it is concluded that the data of the constant amplitude 
load sequences (SOI and S03) are described properly. The spike load sequence (S06) is more 
difficult. In the lower AK regime the crack growth rate is underestimated by the CORPUS model. 

At higher AK Levels the contribution of the primary plastic flow component to the crack growth 
rate is described in a correct way [14] and the value for the yield parameter a NEW in this study can 
be determined. 

Clearly, the addition of the spike load sequence to the constant amplitude sequences helps to 
determine a useful value for c* NEW . Further, it is noted that the more common non-retardation 
models for crack growth prediction can not predict a lower crack growth rate for case S06 
compared to SOI. The contribution of the spikes will slightly increase the crack growth rate in these 
models. From figure 7, however, it can be seen that the crack opening models predict a crack 
growth rate that is one order of magnitude lower (case S06) compared to case SOI. This agrees with 
the test results. 

The values obtained for the material parameters are shown in table 2. Other material properties 
are given in reference [5]. The corresponding values of parameters in the NASGRO crack growth 
law [4] were determined by fitting the same data points. In this case the parameters q and K c are 
thought to represent the static crack growth properties; values for these two parameters are 
calculated directly from m and K ref in the incremental crack growth model discussed here. Further, 
the threshold parameter is kept at the same value, so, in this case the fit parameters are C, n, p and 
"new* * n sequence S06 the contribution of static growth due to the spikes is accounted for in the 
same way as described for the incremental crack growth law. The parameter values obtained are 
shown in table 3. 

Table 2 Material parameter values determined for the incremental crack growth law discussed in 
this study (in MPa, mm) 



7075-T73 

Ti-Al-4V 

2219-T851 

Elber coefficient C 

0.276E-10 

0.294E-10 

0.32527E-10 

Elber exponent n+1 

2.999 

2.656 

2.895 

Threshold exponent p 

1.36 

5.955 

5.948 

Threshold level AK^ 

87.49 

506.28 

249.20 

Static growth exponent m 

17.64 

10.06 

12.99 ! 

i Static growth parameter K,. ef 

1945.2 

4427.2 

1889.4 | 

Yield parameter a NEW 

1.053 

1.197 

1.24 
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7 . 7a Plot of the crack growth equation fitted to 
the measured data points. The material is 
Ti-6AI-4V. The thickness is 2 mm. 
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Fig. 7b Plot of the crack growth equation fitted to 
the measured data points. The material is 
7075-T7351. The thickness is 2 mm. 
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Fig. 7c Plot of the crack growth taw fitted to 

the measured data points. The material is 
2219T851. The thickness is 2 mm. 
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Table 3 Material parameter values determined for the NASGRO modified Forman crack growth 
law (MPa, mm) 



7075-T73 

Ti-Al-4V 

2219-T851 

Elber coefficient C 

0.1745E-9 

0.2621E-10 

0. 135E-10 

Elber exponent n 

2.657 

2.656 

3.399 

Threshold exponent p 

1.150 

0.113 

0.100 

Threshold level AK^ 

87.49 

506.28 

249.20 

Static growth exponent q 

0.208 

0.382 

0.213 

; Critical K c 

2063.30 

5254.92 

2040.89 

Yield parameter a NEW 

1.049 

1.199 

1.358 


Results obtained for the sheet specimens 

Using the NAGRO programme and the values of the material parameters from table 3 the crack 
growth rates and crack size were predicted. The results are collected in figures 8 and 9. Each plot 
represents one test case that includes the experimental result and the result predicted by the 
CORPUS and the STRIP- YIELD models. The first symbol in the identification indicates the type of 
material (Sheet or Plate). The next digits give, the sequence number of the test and the rest 
identifies the material. 


Comparison of dc/dN versus c data (Figs. 8a and 8b) 

The data sets used for the determination of the material parameter values were SO 17075, 
S037075 and S067075. From the corresponding plots it is seen that the "fit" results are reproduced 
properly by the CORPUS and by the STRIP- YIELD modules in NASGRO. This justifies the use of 
the CORPUS crack opening model for determination of the material parameters. The rather strong 
deviations observed for smaller crack sizes in the S05, S06 and S07 spike load tests are probably 
due to the simple definition of the constraint parameters. The same applies to all crack sizes for the 
relatively low spike level in test S14. It appears that too much plastic deformation is predicted by 



crack length c (mm) crack length c (mm) 

Fig. 8a Comparison of measured and predicted crack growth rates, SHEET specimens 
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Fig. 8b Comparison of measured and predicted crack growth rates, SHEET specimens (continued) 
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the STRIP- YIELD model. So, the constraint effect on yielding in secondary flow in tension is 
underestimated. In constant amplitude loading without spikes this effect is absent. The data 
predicted by the STRIP- YIELD model for constant amplitude loading at R = 0.70 are based on a 
cycle by cycle analysis. This guarantees that the contact area at the minimum load is properly 
modelled (this area is of the same order of magnitude as the element size applied in front of the 
crack tip). 

In general, the first cycle of the load sequences shows a different crack growth rate. This is caused 
by the model used to predict growth during primary plastic flow. The deviations disappear if the 
average values for the crack growth rate are plotted in stead of the values per load cycle. It is also 
observed that the STRIP- YIELD results show some transient behaviour at the beginning of the 
curves. The CORPUS model does not account for such effects. The interval between two measured 
data points was to large to see a similar behaviour in the tests. 

The three times repeated single spike sequences, S08 and SI 1, were included in the test program 
to verify the effect of constraint on the retardation region. As this should be judged in a c versus N 
plot these results are discussed in the next section. 


The behaviour of the specimens subjected to the STS load spectrum (SI 3) are also shown. The 
differences are fairly large. The origin of the differences is unknown. For comparison the results 
obtained for a non-retardation model are also presented. 


From the results discussed so far it is concluded that the definition of the constraint factors 
needs improvement. Further, an increased value of the constraint factor for secondary plastic flow 
in tension must be used. This will certainly improve the accuracy of the spike load cases as for 
smaller crack sizes the predicted crack growth rates will increase. For larger crack sizes the effect 
will diminish compared to the contribution of static crack extension to the crack growth rate. In 
view of these observations a judgement on the basis of crack size versus number of cycles would be 
cripple as the deviations observed for the smaller crack sizes will govern the whole picture. For this 
reason this judgement is postponed until the constraint problems are solved in a satisfactory way. 

Comparison of c versus N data (Fig. 9) 


The results obtained for the sheet specimen subjected to three single spikes are given in 


c versus N for si 17075 



Fig. 9 Comparison of measured and predicted crack 
growth curves. SHEET specimen subjected to 
C.A. loading and three single spike loads 
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figure 9. It is seen that the predicted effects of overloads on the crack growth rate is overestimated. 
The 1.5 overload level (SI 1) gives less retardation and, as a result, all curves are more close to the 
case of constant amplitude loading. The spike loads are applied at approximately the same number 
of load cycles as used in the experiment. This implies that the crack sizes (and K factors) are 
smaller and so the plastic zone size and retardation region. As most of the spikes are applied at 
crack sizes smaller than 10 mm (pure plane stress conditions not yet present) it must be concluded 
that under plane strain conditions the constraint factor a is higher than assumed in the models. The 
results obtained for specimen S08 were out of range. These results are left out of the discussion. It 
is noted that predictions based on linear models do not show any retardation. 

Results obtained for the plate materials 

Comparison of dc/dN versus c data (Fig. 10a and 10b) 

The predictions for the plate material are based on the same material parameter values as used 
for the sheet specimens. By nature, the influence of bending stresses induced by clamping the thick 
specimens is more pronounced than for the case of sheet specimens. This may explain part of the 
larger differences in results obtained for the plate specimens. 

For the spike load sequences the conclusion can be drawn that the crack opening load is too 
high. As mentioned previously this indicates that too much plastic deformation occurs in the 
constant amplitude cycles and gives support to the introduction of an increased value for the 
constraint factor for secondary cyclic plastic flow in tension. 

The prediction of the specimen under STS spectrum loading is conservative. Further, the 
deviations between the CORPUS and the STRIPY predictions are relatively large. The source is 
unknown. Again, it is concluded that the definition of the constraint factor a needs reconsideration. 


Comparison of c versus N data (Fig. 11) 

The behaviour observed for the single spikes applied three time on the plate specimens is shown 
in figure 11. The following conclusions can be drawn. The agreement between measured data points 
and predictions is good for the 1.5 overload cases (P05). For the 1.7 spike load case the predicted 
delay is far to large for the case (P05), indicating that the constraint factor a must have a higher 
value to reduce the plastic zone size and the retardation region. 



Fig. 10a Comparison of measured and predicted 
crack growth rates. PLATE specimen 
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Fig. 1 0b Comparison of measured and predicted crack growth rates, PLA TE specimens (continued) 
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Fig. 1 1 Comparison of measured and predicted crack growth curves, PLATE specimens 
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CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 


1 . The empirical CORPUS model and the mechanical STRIP- YIELD model for prediction of the 
crack opening load level were designed and implemented in the NASGRO software to yield the 
NASGRO-STRIPY-93 software. This software was tested successfully. 

2. Material parameter values were determined successfully using two constant amplitude load 
sequences executed respectively at R values of 0.70 and 0.05. In addition a repeated spike load 
sequence was introduced to determine the value of the yield parameter <* NEW ar *d the parameters 
involved in static crack extension during primary plastic flow. For all three materials, Ti-6A1- 
4V (t = 2 mm), 7075-T73(5 1) (t ~ 2 mm) and 2219T851 (t = 4 mm) the material parameters 
were determined. 

3. The model used to describe threshold behaviour under variable amplitude loading is sufficiently 
accurate. In eq. (1) the threshold parameter C 2 is still to be determined. Further research is 
recommended. 

4. The model introduced to describe accelerated crack growth during primary plastic deformation 
in virgin material describes the high crack growth rates observed at high K 1Iiax levels in a 
satisfactory way. The iterative scheme for computation of the transition levels K* appears to be 
functioning properly. 

5. It is important to note that constraint effects primarily depend on the state of stress and, 
therefore, on the loading history, as such a discussion in terms of material parameters is not 
very useful. 

6. The definition and quantification of constraint factors need further improvement. The best way 
to proceed is to collect and evaluate full 3D Finite Element analyses results of the elastic-plastic 
deformation and crack growth problem. These are not yet available in sufficient detail. For the 
time being we need engineering judgement and empirical methods to proceed. 

7. The results obtained for constant amplitude loading (cases SOI, S03 and P01) are covering states 
of stress ranging from plane strain to plane stress conditions. This indicates that in constant 
amplitude loading the crack opening level primarily depends on the yield ratio parameter c* NEW . 
The results obtained for constant amplitude loading interrupted for application of single spikes 
indicate that under plane strain conditions the yield stresses in tension are much higher than 
assumed in this study (the value of appears to be independent of the state of stress). Thus, 
the constraint system depends on the state of stress. In addition the transition from plane strain 
to plane stress must be modelled. This can be done in a way as described in references [9] and 
[ 11 ]. 

8. It is not likely that the constraint effects at different locations in the wake of the crack (resulting 
from residual stresses) are the same as the constraint effects active in yielding in compression in 
front of the (closed) crack. Measurement of the residual stresses on the fracture surface is 
strongly recommended to quantify these constraint effects. 
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9. The results of this verification indicate that the accuracy of the CORPUS model is close to the 
accuracy of the STRIP-Yield model in its current formulation. However, the possibility to 
improve the description of constraint effects gives the STRIP-YIELD model the potential to be 
superior compared to the empirical CORPUS model. Further, compared to the common non- 
retardation models for crack growth prediction a better description of the behaviour observed for 
spike load sequences is demonstrated for both the CORPUS and the STRIP-YIELD models. 
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APPENDIX A The NLR/CORPUS crack opening function 


The crack opening function f(a NEW ,R,S max /<r) 
following way: 



S 

max 


is defined in the 


f - 1 «NEW + "NEW * CF1 






1 

3 

- 1 

J 

'X 

CF1 + _ 
2 

1 

1 

"new 

"new j 


CF2 


f = — CF2 + 
"new 


1 

3 

"NEW j 


CF3 


0 < "new < 1 

1 < a NEW ^ 3 

"NEW > 3 


In the current CORPUS module the functions CF are assumed to be the following form: 


CF1 = (.505 + .18R - . 135R 2 + .81R 3 - .36R 4 ) * Hr 
CF1 = (.505 + . 18R) * Hr 

CF2 = .25 + .06R 4- 1.13R 2 - .44R 3 
CF2 = .25 -F .06R 

CF3 = R 
CF3 = 0 

where: 


Hr 


1 - 0.2(1 - R) 3 


S 


max 


1 3 


l ^ J 


0 <_ R < 1 
R < 0 

0 j< R < 1 
R < 0 

0 <. R < 1 
R < 0 


374 



